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Abstract 

We study the effect of a magnetic field on the thermal relaxation of non-interacting 
small monodomain particle systems particles with a distribution of anisotropy con- 
stants and random easy-axes directions. Numerical calculations of the relaxation 
curves for different distribution widths, and under different magnetic fields H and 
temperatures T, have been performed in the framework of a two-state approxima- 
tion. We show how the obtained data can be analyzed in terms of an modified 
Tln(t/To) scaling from which the field dependence of the mean relaxing energy bar- 
riers can be extracted, a microscopic information which is not easily obtainable by 
other methods. 
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1 Introduction 



The progressive miniaturization of the elementary recording units for magnetic 
recording has put forward the problem of stability of magnetization against 
thermal activation. [1,2] Therefore, the study of the thermal relaxation of 
the magnetization in fine particle systems is of crucial interest to ascertain 
the stability of the recorded information. Although the first studies in this 
field date back to the 40 's [3-5], at present, there is no complete model to 
account for all the experimental findings. One of the points that has not been 
completely clarified is the influence of a magnetic field in the relaxation. 

We want to account for the experimental studies on the relaxation of small- 
particle systems, which essentially measure the acquisition of magnetization 
of an initially demagnetized sample under the application of a magnetic field. 
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[6-10] In this kind of experiments, the field modifies the energy barriers of 
the system that are responsible for the time variation of the magnetization, 
as well as the final state of equilibrium towards which the system relaxes. 
Here, we will demonstrate that the previously proposed Tln(t/r ) scaling of 
the relaxation data at zero field [6,11-16] can also be extended to systems 
relaxing in an applied magnetic field. We will show how information about 
the microsopic energy barriers responsible for the relaxation can be obtained 
from this procedure. 

The paper is organized as follows. In Sec. 2, we present the basic features 
of model showing how the distribution of energy barriers of the system is 
influenced by the application of a magnetic field with the help of the concept of 
effective energy barrier distribution. In Sec. 3, we derive the equation governing 
the time dependence of the magnetization from a master rate equation for a 
Two-State Approximation (TSA). The results of numerical calculations based 
on the above mentioned equation are presented in Sec. 4. There, we present the 
Tln(t/ro) scaling of relaxation curves at a given magnetic field, discussing its 
range of validity. We also study the possibility of a scaling at different fields 
and fixed temperature, and its applications. Finally in Sec. 5 we resume the 
main conclusions of the article. 



2 Model 



We consider an ensemble of randomly oriented, noninteracting single-domain 
ferromagnetic particles of volume V and magnetic moment M = M s Vm with 
uniaxial anisotropy. To take into account the spread of particle volumes in 
real samples, we will assume that the particles anisotropy constants K are 
distributed according to a logatithmic-normal distribution 

/(if) = * e -lnHK/Ko)/2v* (1) 



The energy of a particle is determined by the orientation of M with respect 
to the external magnetic field H and to the easy-axis direction n. Using the 
angular coordinates defined in Fig. 1, it can be written as 

E = = — cos 2 (#) — 2h cos(# — ip) . (2) 



where we have defined the reduced field h = H/H c and H c = 2K/M S as 
the critical field for an aligned particle. We have concentrated on the two 
dimensional case (M lying in the plane formed by H and n) since the energy 
maxima and minima can be calculated analytically only in this case. In Fig. 
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Fig. 1. Energy function E(6, ip) as a function of the angle between the magnetization 
vector m and the magnetic field h, for m in the plane formed by the easy-axis and 
the magnetic field (ip = 0), as given by Eq. (2). The plot is for a particle whose 
easy-axis n forms an angle ip = 30° with h, and H = 0.3. We have used the 
following notation to designate the extrema of the energy: 9^ ia and # max refer to 
the extrema closer to the field direction while 6?l- and 0!L V refer to those farther 

mm nidx 

from the direction of the field. The four possible energy barriers between them are 
E^ J = i?(# max ) — E(9 :) m - m ). Inset: Schematic representation of the quantities involved 
in the definition of the system. The easy-axis of the particles n are in the x-z plane 
forming an angle ip with the magnetic field H, which points along the z axis. 6 and 
ip are the spherical angles of the magnetization vector M. 
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1, we show the variation of the energy with 9 for a typical case, defining in 
the same figure the notation for the energy barriers and extrema. 



The magnetic field modifies the energy barriers of the system depending on 
the particle orientation and anisotropy value, and, consequently, changes the 
original energy barrier distribution [17,18]. In order to understand the qual- 
itative change of f e s with h, we have numerically calculated f^E^) for a 
system of randomly oriented particles with f(K), for different widths a and 
K — 1, and several values of the magnetic field h. The calculation has been 
performed by making energy barriers histograms for a collection of 10 000 
particles. The results are given in Fig. 2. In all the cases, we observe the pro- 
gressive splitting of the original distribution f(E®) in two subdistributions of 
high and low barriers as h increases from zero. The field tends to make deeper 
one of the minima, therefore increasing the two energy barriers for rotation 
of M out of the field direction, while the other two are reduced. In this way, 
the global effect of h is a splitting of f(E^) towards lower and higher values of 
E\>. As h attains the critical value h c for the particles with smaller K, a peak 
of zero or almost zero energy barriers starts to appear (see for example the 
curves for H = 0.5 in the case a = 0.5); while most of the non-zero barriers 
are distributed according to a distribution identical to f(E®), but centered 
at higher energies. The higher the width of the distribution a, the lower the 
h at which the lowest energy barriers start to be destroyed by the field. The 
random orientation of the anisotropy axes smears out the distributions and 
makes the minima less pronounced than for a system with aligned particles 
due to the spread in particle orientations, this effect being more noticeable at 
high h. In Sec. 4, we will discuss how these results affect the time dependence 
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Fig. 2. Effective energy barrier distributions for particles with random orientations 
of anisotropy axes and a lognormal distribution of anisotropy constants of a = 0.2 
(continuous lines), a = 0.5 (dashed lines), for values of H as indicated in the figures. 
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of magnetization in relaxation experiments. 



3 Relaxation curves in the presence of a magnetic field 



In order to derive an equation for the time dependence of the magnetization 
we propose a simple dynamical model that can be solved analytically in terms 
of intrinsic parameters. Since usually, in experiments performed in real sam- 
ples the thermal energy ksT is smaller than the anisotropy energy KV, we 
will concentrate on the high energy barrier (low T) regime. In this case, it is 
justified to consider only transitions between the two energy minima of each 
particle, and to truncate the continuum of possible states for the magnetiza- 
tion to the two energy minima [19,20] (TS A). The master equation for Pj, the 
probability to find the magnetization of the particle in the minimum i, can be 
written as [21] 

f-EEW'i-*}, (3) 



(k) 

where stands for the transition rate for a jump from the state % to the state 
j separated by the maximum k (see Fig. 1). Our choice for these quantities 
takes into account the energy barrier separating the initial minimum % from 
the final state j 

(*) J_ 1_ E k if) 
w l3 (fc) e , W 



where r 1 is the attempt frequency and (3 = 1/ksT. 

Taking into account the normalization condition Pi + P 2 — 1, one can easily 
solve Eq. (3) for Pi and P 2 as a function of time 



1 _ pPsp-t/r 

w 1 + eP £ K J 

The time-dependence of the system is thus characterized by an exponential 
function with a single relaxation time r that takes into account all possible 
probability fluxes 

r- 1 = To" 1 (e-^ 2 + e"^ 2 ) (l + e^) , (6) 
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where e = E^ in — Ef nin is a measure of the asymmetry of the energy function. 

The time dependence of the magnetization of the particle is then finally given 
by: 

m{t-K^) = cos[^ in (^)]Pi(t) +cos[e 2 min (4>)]P 2 (t) 

= rh TS (K,^) + [m - m TS (K, ^e'^^. (7) 

In this equation, t?its(K, ip) is the equilibrium magnetization in the TSA, and 
m is the initial magnetization. If we have an ensemble of randomly oriented 
particles and a distribution of anisotropy constants f(K), then the relaxation 
law of the magnetization is given by 

OO 7T 

m(t) = J dKf(K) J dip m(t;K,iP) . (8) 
o o 

This will be the starting point for all the subsequent numerical calculations 
of the relaxation curves and magnetic viscosity. 



4 Numerical calculations 

4-1 Relaxation curves: Tln(t/r ) scaling and normalization factors 

In this section, we present the results of numerical calculations of the magne- 
tization decay based on Eq. (8). For the sake of simplicity, we have assumed 
zero initial magnetization m = 0, so particles have initially their magnetic 
moments at random and evolve towards the equilibrium state m cq . In the fol- 
lowing, we will use dimensionless reduced variables for temperature and time, 
defined as T/T and t/r , with T = E^/k-Q and E the value of the energy at 
which f(K) is centered. In Fig. 3, we show the results of the numerical calcu- 
lations for a system with a = 0.5 and H = 0.5 and temperatures ranging from 
0.02 to 0.2. In the left hand panel, we present the original relaxations nor- 
malized to the equilibrium magnetization value that has be computed in the 
TSA. Normalization is essential in order to compare relaxations at different 
temperatures [14], especially at low fields where the temperature dependence 
of the equilibrium magnetization is more pronounced. 

Our next goal is to investigate the possibility of scaling relaxation curves at 
different T in a given magnetic field with the scaling variable T ln(t/ro), in the 
spirit of our previous works [11-14]. Let us briefly recall the basic points of 
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Fig. 3. (a) Relaxation curves for an ensemble of particles with randomly oriented 
anisotropy axes and the distribution f(K) of Eq. 1 of width a = 0.5 and Kq = 1 
calculated by numerical integration of Eq. (8). Reduced temperatures T/Tq, starting 
from the lowermost curve, range from 0.01 to 0.1 with 0.01 increments, and from 0.1 
to 0.2 with 0.02 increments. The applied field is H = 0.5. The original relaxation 
curves have been normalized to the equilibrium magnetization mT,s(T). (b) Relax- 
ation curves of panel (a) plotted as a function of the scaling variable Tm(t/ro) to 
obtain the master relaxation curve. 

the method. The time dependence of the remanent magnetization of a system 
with a distribution of energy barriers f{E) can be written as 

oo oo 

M(t) = M J dEf(E) e-* /r(B) ~ M J dEf(E) , (9) 

E c (t) 



with a characteristic relaxation time given by the Arrhenius-Neel theory as 
t(E) = r exp E / T . Then, the relaxation law for the magnetization will be a 
function of the parameter E c (t) = Tln(£/ro) if the Boltzmann exponential 
factor in (9) can be approximated by a step function centered at E c (t). Then, 
as was shown in [11], the requirement for the validity of this approximation 
is that the characteristic width of f(E) is of the order of the width of the 
step-like function e~ f ^ E \ which can be calculated as A = eT (e is the Napier 
number). This implies low enough T or energy distributions spread over a 
wide range of energies (large a). 

In Fig. 3b, we show the original relaxation curves of Fig. 3a as a function 
of the scaling variable Tln(t/r ). The collapse of the curves into a single 
master curve indicates that scaling is accomplished also in the presence of 
a magnetic field. Only at short times, the higher T curves depart from the 
master curve because in this T range the scaling requirements are no longer 
accomplished. According to the requirements for Tln(t/r ) scaling presented 



7 



in the preceding paragraph, in the presence of a magnetic field, the validity of 
the scaling hypotehsis at a given T depends on the width of the effective energy 
barrier distribution (see Fig. 2). From simulations of master curves obtained 
at different H (not presented in Fig. 3), we observe that, the higher the field, 
the better the scaling of the curves is in the long time region and the worse at 
short times. This observation holds independently of the value of a, indicating 
that it is a consequence of the application of a magnetic field. This can be 
understood with the help of the effective energy barrier distribution introduced 
in Sec. 2. As was shown in Fig. 2, H widens f e s{E) and shifts the lowest 
energy barriers towards the origin, giving rise to a subdistribution of almost 
zero energy barriers that narrows with increasing H, and, consequently, the 
requirements for Tln(t/ro) scaling are worse fulfilled at small T ln(t/ro) values. 
On the contrary, as itwas shown in section 2, H broadens the high energy tail 
of energy barriers that contribute to the relaxation, f(E^), improving the 
scaling requirements at large Tln(t/r ) values. 

4-2 Scaling of relaxation curves at different magnetic fields 

Another interesting point is the possibility of finding an appropriate scaling 
variable to scale relaxation curves at different fields for a given T, in a way sim- 
ilar to the case of a fixed field and different temperatures, in which Tln(t/r ) 
is the appropriate scaling variable. In a first attempt, we will study the effect 
of h on a system with random anisotropy axes and distribution of anisotropics 
f(K) of width a, Eq. (1). Low temperature relaxation rates, obtained by per- 
forming the logarithmic time derivative of M(t) 

s «=liH^- T »- < 10 > 

which is the coeficient of magnetic viscosity S(t) are presented in Fig. 4, for 
o = 0.2, 0.5. From this expression, we see that when S is evaluated by taking 
the derivative of the master relaxation curve with respect to T\n(t/r ), the 
resulting curve is proportional to the effective energy barrier distribution. 
First, let us notice that the qualitative shape of the viscosity curves is not 
distorted by H. The magnetic field simply shifts the position of the maxima 
towards lower values of Tln(t/ro) and narrows the width of the peaks, these 
effects being similar for both studied o. The energy corresponding to the 
maximum relaxation rate (E™ ax ) decreases with increasing h, following the 
decrease of the smallest energy barriers (see Fig. 5), with an almost linear 
dependence on h. As can be clearly seen in Fig. 5, the dependence of E max on 
h (in symbols) follows that of the lowest energy barriers for particles oriented 
at ip = 7r/4 and with K = K max (dashed lines), which have been evaluated 
numerically from Eq. 2. Let us also notice that .Smax becomes almost constant 
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below ho (the value at which the field starts to destroy the lowest energy 
barriers) and low enough T (see the position of the maxima in the curves of 
Fig. 4). indicates that the relaxation curves at different h and fixed T may 
be brought to a single curve by shifting them along the Tln(t/r ) axis in 
accordance to the E max variation. The resulting curves are displayed in Fig. 6 
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Fig. 4. Low temperature (T = 0.05) viscosity curves for a system of particles with 
random orientations and lognormal distribution of anisotropies with (a) a = 0.2 and 
(b) a = 0.5. The curves have been normalized to the equilibrium magnetization and 
correspond to magnetic fields H = 0.1 to 1.0 in 0.1 steps starting from the right. 
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Fig. 5. Field dependence of the energy corresponding to the maximum relaxation 
rate, T ln(t max / Vo) , as derived from the viscosity curves in Fig. 4 for T = 0.05 
(circles). The dashed line shows the h dependence of the lowest energy barriers for 
the particles oriented at ip = 7r/4 and K equal to the most probable value. 



9 







0.5 



1.0 



1.5 



2.0 




Fig. 6. Normalized relaxation curves as a function of the scaling variable 
Tln[(t/t max (h)] for T = 0.05) obtained from Fig. 4 by shifting the curves in the 
horizontal axis with the position of the maximum relaxation rate of Fig. 5. For 
a = 0.2 the field are H = 0.1, 0.2, 0.3, 0.4, 0.5, and for a = 0.5, H = 0.1, 0.2, 0.3, 0.4 
(starting from the uppermost curve). 

for a = 0.2, 0.5. They are the equivalent of the master curves of Fig. 3 for a 
fixed h and different T. Now the appropriate scaling variable is 



which generalizes the scaling at fixed T (another approach to a similar scaling 
law may be found in [22]). This new scaling is valid for fields lower than ho, 
the field at which the lowest barriers start to be destroyed and above which 
the relaxation becomes dominated by almost zero energy barriers. Thus, as 
already discussed in the previous paragraphs, the wider a, the smaller the h 
range for the validity of field scaling. 



5 Conclusions 

We have proposed a model for the relaxation of small particles systems under 
a magnetic field which can be solved numerically from analytical expressions 
and which allows to study the effect of the magnetic field on the energy barrier 
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T\n[t/t max (h)} , 



(11) 
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distribution. In particular, we have shown that the original energy barrier dis- 
tribution f(E-b) is split into two subdistributions which evolve towards higher 
and lower energy values, respectively, as h increases. 

For fields smaller than the critical values for the smallest barriers, the relax- 
ation curves at different h and fixed T can be collapsed into a single curve, in 
a similar way than Tln(t/ro) scaling for curves at fixed h. Whereas the latter 
allows to extract the barrier distribution by differentiation of the master curve 
[14], the shifts in the T\n(t/r ) axis necessary to produce field scaling, give 
the field dependence of the mean relaxing barriers, a microscopic information 
which cannot easily be inferred from other methods [23]. 
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